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X Abstract 

A search for the single production of doubly-charged Higgs bosons H ±i: in ep collisions 
is presented. The signal is searched for via the Higgs decays into a high mass pair of 
same charge leptons, one of them being an electron. The analysis uses up to 118 pb _1 of 
ep data collected by the HI experiment at HERA. No evidence for doubly-charged Higgs 
production is observed and mass dependent upper limits are derived on the Yukawa cou- 
plings hd of the Higgs boson to an electron-lepton pair. Assuming that the doubly-charged 
Higgs only decays into an electron and a muon via a coupling of electromagnetic strength 
h e/1 = y/47ra em = 0.3, a lower limit of 141 GeV on the H ±ziz mass is obtained at the 95% 
confidence level. For a doubly-charged Higgs decaying only into an electron and a tau and 
a coupling h eT = 0.3, masses below 112 GeV are ruled out. 
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1 Introduction 

Doubly-charged Higgs bosons (H ±=t ) appear when the Higgs sector of the Standard Model 
(SM) is extended by one or more triplet(s) with non-zero hypercharge [1-3]. Examples are pro- 
vided by some Left-Right Symmetric models [4], or their supersymmetric extensions, which 
are of particular interest since they provide a mechanism to generate small non-zero neutrino 
masses. Such models can lead to a doubly-charged Higgs boson light enough [5] to be pro- 
duced at the existing colliders. The Higgs triplet(s) may be coupled to matter fields via Yukawa 
couplings which are generally not related to the fermion masses. A non-vanishing coupling 
of a doubly-charged Higgs to an electron would allow its single production in ep collisions at 
HERA. This possibility is investigated in this paper with a search for doubly-charged Higgs 
bosons decaying into a high mass pair of same charge leptons, one of them being an electron. 

An analysis of multi-electron events was already presented by the HI collaboration [6]. Six 
events were observed with a di-electron mass above 100 GeV, a domain in which the Standard 
Model prediction is low. In the present paper the compatibility of these events with the hy- 
pothesis of a doubly-charged Higgs coupling to ee is addressed and a further search for a H ±± 
boson coupling to e/i and er is performed. The analysis is based on ep data collected by the HI 
experiment between 1994 and 2000, which amount to a luminosity of up to 118 pb _1 . 



2 Phenomenology 

At tree level, doubly-charged Higgs bosons couple only to charged leptons and to other Higgs 
and gauge bosons. Couplings to quark pairs are forbidden by charge conservation. The cou- 
plings of a doubly-charged Higgs to charged leptons can be generically described by the La- 
grangian: 

C = Y J h^Hi%l l c P L , R l ] +h.c. , (1) 

where / are the charged lepton fields, l c denote the charge conjugate fields, i,j are generation 
indices, and Pl,r = (1 =F 7s)/2. The Higgs fields H~^\ coupling to left-handed or right-handed 
leptons correspond to different particles and not all models predict their simultaneous existence. 
The Yukawa couplings \\ . — h t .\. are free parameters of the model. 

The phenomenology of doubly-charged Higgs production at HERA was first discussed in [7]. 
For a non-vanishing coupling hj the single production of a doubly-charged Higgs boson is 
possible at HERA in ey* interactions via the diagrams shown in Fig. 1, where a photon is 
radiated off the proton or one of its constituent quarks. The proton may remain intact or be 
broken during this interaction, leading to an elastic or inelastic reaction, respectively. With 
longitudinally unpolarised lepton beams, as were delivered by HERA until 2000, the H ±± 
production cross section does not depend on whether the Higgs couples to left-handed or right- 
handed leptons. Hence a generic case is considered here of a doubly-charged Higgs boson 
which couples to either left-handed or right-handed leptons and the L, R indices are dropped in 
the following. 
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Figure 1 : Diagrams for the single production of a doubly-charged Higgs boson in e + p collisions 
at HERA via the h d coupling. The hadronic final state is denoted by p (X) in the elastic 
(inelastic) case, where the initial proton remains intact (dissociates). The contribution of Z 
exchange can be safely neglected. 

Within the mass range considered in this analysis, it is assumed that decays of the H ±± into 
gauge bosons and other Higgs particles are not allowed kinematically such that the doubly- 
charged Higgs only decays via its Yukawa couplings into a lepton pair. 

Indirect upper bounds on the Yukawa couplings of a doubly-charged Higgs to leptons are 
reviewed in [8]. The coupling h ee of a doubly-charged Higgs to an electron pair is con- 
strained by the contribution of virtual H ±± exchange to Bhabha scattering in e + e~ collisions. 
A recent OPAL analysis [9] sets the constraint h ee < 0.14 for a doubly-charged Higgs mass 
M H = 100 GeV. From low energy e + e~ data, coupling values of 0(0.1) are allowed for h eil 
and h eT for a Higgs mass of 100 GeV [10]. Taking these indirect constraints into account, 
the production of a doubly-charged Higgs mediated by h ee , h eil or h eT might be observable at 
HERA. The Higgs signal would manifest itself as a peak in the invariant mass distribution of 
same charge ee, e/i or er leptons, respectively. For the range of masses and couplings probed 
in this analysis, the Higgs decay length is vanishingly small but its width remains negligible 
compared to the experimental resolution on the mass of the lepton pair. 



3 Simulation of the Signal and Standard Model Backgrounds 

The calculation of the cross section for doubly-charged Higgs production, as well as the simula- 
tion of signal events, relies on a dedicated Monte Carlo program developed for this analysis. The 
differential cross sections are integrated using the VEGAS package [11]. Different approaches 
are followed depending on the photon virtuality Q 2 and on the mass W of the hadronic final 
state: 

• in the inelastic region (W > m p + m^, with the proton mass m p and the pion mass 
m^) and when the photon virtuality is large (Q 2 > 4 GeV 2 ), the interaction involves 
a quark inside the proton. The squared amplitude of the process e ± q — ► e T H ±± q is 
evaluated using the CompHEP package [12, 13]. The parton densities in the proton are 
taken from the CTEQ4L [14] parameterisation and are evaluated at the scale \[Q 2 . The 
parton shower approach [15] based on the DGLAP [16] evolution equations is applied to 



simulate QCD corrections in the initial and final states, and the hadronisation is performed 
using PYTHIA 6.1 [15]. 

for the elastic region (W = m p ) and the inelastic region at low Q 2 (W > m p + m^, 
Q 2 < 4 GeV 2 ), the squared amplitude is calculated using the FORM program [17]. The 
hadronic tensor is parameterised in terms of the usual electromagnetic structure functions 
F\(x, Q 2 ) and F 2 (x, Q 2 ) of the proton, where x = Q 2 /(W 2 + Q 2 — rn 2 ). For the elastic 
process these structure functions are expressed in terms of the electric and magnetic form 
factors of the proton. For the low Q 2 inelastic region they are taken from analytical 
parameterisations [18]. The simulation of the hadronic final state for low Q 2 inelastic 
events is performed via an interface to the SOPHIA program [19]. 



For a Yukawa coupling h ee or h efl of electromagnetic strength (h = ^/4ira em = 0.3) the total 
cross section amounts to 0.39 pb (0.04 pb) for a Higgs mass of 100 GeV (150 GeV). The low 
Q 2 (high Q 2 ) inelastic contribution is found to be ~ 30% (~ 20%) of the total cross section in 
the mass range 80 — 150 GeV. The cross section for producing a doubly-charged Higgs via a 
coupling h eT is lower by about 40% due to the non-negligible mass of the r lepton produced in 
association with the Higgs. 

The theoretical uncertainty on the cross sections obtained is taken to be 4% in the mass range 
considered. This is derived from an assessed uncertainty of 2% on the proton form factors [20] 
and from the uncertainty on the scale at which the parton densities for the inelastic contribution 
are evaluated. The latter uncertainty is estimated from the variation of the computed cross 
section as this scale is changed from a/Q^/2 to 2^/Q 2 . 

Separate signal event samples corresponding to the production and decay of a doubly-charged 
Higgs via a coupling h ee , h etl and h eT have been produced for Higgs masses ranging between 
80 and 150 GeV, in steps of 10 GeV. 

Di-electron production, which proceeds mainly via two-photon interactions, constitutes an ir- 
reducible background for ee final states. The production of muon or tau pairs constitutes a 
background for the efi and er analyses when the scattered electron is detected. The Cabibbo- 
Parisi process ee — > 7, Z — ► 11, in which the incoming electron interacts with an electron 
emitted from a photon radiated from the proton, contributes at high transverse momentum only. 
The Drell-Yan process was calculated in [21] and found to be negligible. All these processes 
are simulated using the GRAPE Monte Carlo generator [22], which also takes into account 
contributions from Bremsstrahlung with subsequent photon conversion into a lepton pair and 
electroweak contributions. 

Experimental backgrounds come dominantly from Neutral Current Deep Inelastic Scattering 
(NC DIS) where a jet is misidentified as an electron, a muon or a tau. Compton scattering is 
also a source of background for ee final states when the photon is misidentified as an electron. 
These processes are simulated with the DJANGO [23] and WABGEN [24] generators. 

All generated events are passed through the full simulation of the HI apparatus and are recon- 
structed using the same program chain as for the data. 



4 The HI Detector 



A detailed description of the HI experiment can be found in [25]. Only the HI detector 
components relevant to the present analysis are briefly described here. Jets and electrons 
are measured with the Liquid Argon (LAr) calorimeter [26], which covers the polar angle 1 
range 4° < 9 < 154°. Electromagnetic shower energies are measured with a precision of 
a(E)/E = 12%/^/E/GeV®l% and hadronic energies with a(E)/E = 50%/y/E/GeV®2%, 
as determined in test beams [27]. In the backward region a lead/scintillating-fibre 2 (SpaCal) 
calorimeter [28] covers the range 155° < 9 < 178°. The central (20° < 9 < 160°) and for- 
ward (7° < 9 < 25°) tracking detectors are used to measure charged particle trajectories, to 
reconstruct the interaction vertex and to supplement the measurement of the hadronic energy. 
The LAr and inner tracking detectors are enclosed in a super-conducting magnetic coil with 
a strength of 1.15 T The return yoke of the coil is the outermost part of the detector and is 
equipped with streamer tubes forming the central muon detector (4° < 9 < 171°). In the for- 
ward region of the detector (3° < 9 < 17°) a set of drift chamber layers (the forward muon 
system) detects muons and, together with an iron toroidal magnet, allows a momentum mea- 
surement. The luminosity measurement is based on the Bethe-Heitler process ep — » ep'j, where 
the photon is detected in a calorimeter located downstream of the interaction point. 



5 Data Analysis 

The analyses of ee and efi final states use the full e ± p data set recorded in the period 1994- 
2000, which corresponds to an integrated luminosity of 118 pb _1 . The analysis of er final states 
makes use of the e + p data collected in the years 1996-1997 and 1999-2000, which amount to 
a luminosity of 88pb _1 . The HERA collider was operated at a centre-of-mass energy y/s of 
300 GeV in 1994-1997 and of 318 GeV in 1998-2000. 

Events are first selected by requiring that the longitudinal position of the vertex be within 35 cm 
around the nominal interaction point. In addition topological filters and timing vetoes are ap- 
plied to remove background events induced by cosmic showers and other non-ep sources. The 
main triggers for the events are provided by the LAr calorimeter and the muon system. 

5.1 Lepton Identification 

An electron 3 candidate is identified by the presence of a compact and isolated electromagnetic 
energy deposit above 5 GeV in the LAr or SpaCal calorimeter. The energy of the electron 
candidate is measured from the calorimetric information. In the angular range 20° < 9 < 150° 
the electron identification is complemented by tracking conditions, in which case the direction 

The origin of the HI coordinate system is the nominal ep interaction point, with the direction of the proton 
beam defining the positive z-axis (forward region). The transverse momenta are measured in the xy plane. The 
pseudorapidity r\ is related to the polar angle 9 by r\ = — In tan(6>/2). 

2 Before 1995 a lead-scintillator calorimeter was used. 

3 Unless otherwise stated, the term "electron" is used in this paper to generically describe electrons or positrons. 



of the electron candidate is given by that of the associated track. Electron candidates in the 
forward region, 5° < 9 < 20°, are required to have an energy above 10 GeV. 

A muon candidate is identified by associating an isolated track in the forward muon system 
or in the inner tracking system with a track segment or an energy deposit in the instrumented 
iron. The muon momentum is measured from the track curvature in the toroidal or solenoidal 
magnetic field, respectively. 

Tau leptons are preselected as described in [29] by requiring a track with transverse momentum 
above 5 GeV measured in the inner tracking detector. The leptonic tau decays r — ► evv and 
r — ► [ivv are reconstructed by matching the selected track to an identified electron or muon. 
Tracks that are not identified as electrons or muons are attributed to hadronic tau decays if at 
least 40% of the track momentum is reconstructed in the LAr calorimeter as matched clustered 
energy. In that case it is moreover required that the track belong to a narrow jet: no other 
track should be reconstructed within 0.15 < R < 1.5 around the track direction, where R = 
y/ 'Arj 2 + Aip 2 with Arj and Aip being the distances in pseudorapidity and azimuthal angle, 
respectively. The transverse momentum and the direction of the r candidate are approximated 
by those of the associated track. 

5.2 Analysis of the H — > ee Decay 

This analysis is based on the published HI measurement of multi-electron production [6]. The 
event selection requires at least two central (20° < 9 e < 150°) electron candidates, one of them 
with a transverse momentum P^ 1 > 10 GeV (ensuring a trigger efficiency close to 100% [30]) 
and the other one with Pj? > 5 GeV. After this preselection, 125 events are observed, in good 
agreement with the SM expectation of 137.4± 10.7. In each event, the two highest P T electrons, 
one of those being possibly outside the central region, are assigned to the Higgs candidate. The 
distribution of their invariant mass M ee is shown in Fig. 2a. At low mass a good agreement is 
observed between data and the SM expectation which is largely dominated by 77 contributions. 
Six events are observed at M ee > 100 GeV, compared to the SM expectation of 0.53 ± 0.08. 

Further selection criteria are then applied, which are designed to maximise the sensitivity of 
the analysis to a possible H ±± signal. The charge measurement of the two leptons assigned 
to the Higgs candidate is exploited. In e + p (e~p) collisions, where H ++ (H ) bosons could 
be produced, events in which one of the two leptons is reliably assigned a negative (positive) 
charge are rejected. The charge assignment requires that the curvature k of the track associated 
with the lepton be measured with an error 5k satisfying | k/8k \ > 2. The precise calorimetric 
measurement of the electron transverse momenta is further exploited by applying an additional 
M ee dependent cut on the sum of the transverse momenta of the two electrons assigned to the 
Higgs candidate. The lower bound is optimised to keep 95% of the signal and varies between 
45 GeV and 120 GeV. This cut suppresses events coming from 77 processes. The efficiency for 
selecting signal events varies from 50% for a H ±± mass of 80 GeV to 35% for a H ±± mass of 
150 GeV. In this mass range the resolution on M ee varies between 2.5 GeV and 5 GeV. 

After these requirements, 3 events are observed at M ee > 65 GeV, in agreement with the SM 
expectation of 2. 45±0. 11 events. Amongst the six events 4 atM ee > 100 GeV, only one satisfies 
the final selection criteria. 



4 Out of these, three do not fullfill the M ee dependent Pt cut, and two do not satisfy the charge requirement. 
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Figure 2: Distribution of (a) the invariant mass M ee of the two highest Pp electrons for multi- 
electron events, (b) the electron-muon invariant mass M ejU , and (c) the electron-tau candidate 
invariant mass M eT . The data (symbols) are compared with the Standard Model expectation 
(histogram). The distributions are shown at the preselection level (see text). 

5.3 Analysis of the H — ► e/x Decay 



Events having one electron and one muon with minimal transverse momenta of Py > 10 GeV 
and P^ > 5 GeV are selected. The polar angle of electron candidates is restricted to 20° < 8 e < 
140° to reduce the large background arising from NC DIS events. The 8 range for muon can- 
didates extends towards low angles, 10° < 8^ < 140°, which increases the efficiency for high 
H ±± masses. The minimum transverse momentum required for electron candidates ensures a 
trigger efficiency close to 100% for these events. After this preselection, 35 data events are 
observed compared to a SM expectation of 29.6 ± 3.4. The distribution of the invariant mass of 
the electron and the muon M eM is shown in Fig. 2b. A good agreement is observed between the 
data and the SM expectation, which is dominated by 77 contributions. 

For the final selection of H — > e/x candidates the charge of the e and /x is exploited using the 
same criteria as used in section 5.2. The efficiency for selecting signal events varies from 55% 



to 40% for a H ±± mass between 80 GeV and 150 GeV. The resolution on M eil varies between 
3 GeV and 8 GeV. For M eM > 65 GeV one event is observed while 4.17 ± 0.44 events are 



expected from the SM. 



5.4 Analysis of the H — > er Decay 

The search for a H ++ boson decaying into er is performed in three final states, depending on 
whether the r decays into an electron, a muon or hadronically (h). Details of this analysis can 
be found in [29]. Events are selected which contain either two electrons (ee), or an electron 
and a muon (e/j,), or an electron and a hadronic r candidate (eh) as defined in section 5.1. The 
two leptons, or the electron and the hadronic-r candidate, should have a transverse momentum 
above 5 GeV, be in the angular range 20° < 9 < 140°, and be separated from each other by 
R > 2.5 in pseudorapidity-azimuth. One of them must have a transverse momentum above 
10 GeV, which ensures a trigger efficiency above 95% in all three classes. For events in the e/x 
class the polar angle of the electron candidate is required to be below 120°. 

A significant amount of missing transverse and longitudinal momentum is expected due to the 
neutrinos produced in the r decays. Events in the ee class are required to have a missing trans- 
verse momentum p™ ss > 8 GeV. For the eh class, which suffers from a large NC DIS back- 
ground, it is required that p™ lss > 11 GeV, that the energy deposited in the SpaCal calorimeter 
be below 5 GeV, and that the variable ^ i E % — P l z , where the sum runs over all visible particles, 
be smaller than 49 GeV. For fully contained events ^ E % — P\ is expected to peak at twice the 
lepton beam energy E = 27.5 GeV, i.e. 55 GeV, while signal events are concentrated at lower 
values due to the non observed neutrinos. In total 6 events are preselected, in agreement with 
the SM prediction of 7.8 ± 1.5. 

In each class, the er invariant mass M eT is reconstructed by imposing longitudinal momentum 
and energy conservation, and by minimising the total momentum imbalance in the transverse 
plane. Tau leptons are assumed to decay with a vanishing opening angle. This method yields a 
resolution of about 4 GeV on the mass M eT . Figure 2c shows the er invariant mass distribution 
of the selected events together with the SM expectation. 

For the final selection, events are rejected if the track associated with one of the Higgs decay 
product candidates is reliably assigned a negative charge, opposite to that of the incoming lepton 
beam. The signal efficiencies depend only weakly on M H . The fractions of simulated H — ► er 
events which are reconstructed in the various classes are given in table 1, for an example mass 
of Mh = 100 GeV. The total efficiency on the signal amounts to about 25%. 

The final event yields are also shown in table 1. Only one event (in the eh class) satisfies the 
final criteria, while 2.1 ± 0.5 events are expected. 



5.5 Systematic Uncertainties 

The systematic uncertainties attributed to the Monte Carlo predictions for the ee analysis are 
detailed in [6]. The dominant systematic uncertainty is due to the electron-track association 
efficiency, which is 90% on average with an uncertainty increasing with decreasing polar angle 
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Event 
class 


H ++ — ► e + r + final selection 


N obs 


Nbckg 


Signal fraction 


e\x 





0.27± 0.02 


6% 


eh 


1 


1.66 ±0.48 


12% 


ee 





0.14±0.04 


7% 


total 


1 


2.07± 0.54 


25% 



Table 1: Number of observed (N b s ) and expected (Nbckg) events in each event class which 
satisfy all criteria to select H ++ — ► e + r + candidates with a mass M er > 65 GeV. The last 
column shows the fractions of the H — ► er Monte Carlo events which are reconstructed in the 
various classes, for a mass of 100 GeV. 



from 3% to 15%. Systematic errors due to the uncertainty on the electromagnetic energy scale 
(known at the 0.7% to 3% level in the central and forward regions of the LAr calorimeter, 
respectively) and on the trigger efficiency (3%) are also taken into account. 

For the e/i analysis, the dominant additional systematic uncertainty is due to the muon identi- 
fication efficiency known within 6% [31]. The uncertainty due to the reconstruction efficiency 
of the central tracking detector for central muons contributes an additional 3%. The muon mo- 
mentum scale is known within 5%, and the trigger efficiency for eji final states is known within 

3%. 

The same systematic uncertainties affect the SM expectations in the ee and e/i classes of the er 
analysis. The uncertainty of the hadronic energy scale in the LAr calorimeter (4%) constitutes 
another source of uncertainty due to the cuts applied on the p™' iss and £\ E — P z variables. For 
the eh event class the dominant uncertainties on the SM expectation, coming mainly from NC 
DIS processes, are due to the uncertainty of 3% of the track efficiency, to that of the hadronic 
energy scale, and to that of the hadronisation model. 

The luminosity measurement leads to a normalisation uncertainty of 1.5%. 

For both the expected signal and the predicted background, the systematic uncertainties result- 
ing from the sources listed above are added in quadrature. 



6 Interpretation 

With the final Higgs selection no significant excess over the SM expectation is observed. Upper 
limits on the H ±± production cross section times the branching ratio for the H ±± to decay into 
one of the analysed final states are derived as a function of the H ±± mass and are shown in 
Fig. 3a. The limits are presented at the 95% confidence level and are obtained using a modi- 
fied frequentist approach [32]. Statistical uncertainties, as well as the influence of the various 
systematic uncertainties on both the shape and the normalisation of the mass distributions for 
signal and background events, are taken into account. The best sensitivity is obtained for a H ±± 
produced and decaying via h efl , with upper limits around 0.05 pb. 
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Figure 3: (a) Upper limits at the 95% confidence level on the H ±± production cross section 
times the branching ratio for the H ±± to decay into ee, ejj, or er, as a function of the Higgs 
mass, (b)-(d) Upper limits on the coupling h e i assuming that the H ±± couples only (b) to ee, 
(c) to en or (d) to er. Regions above the curves are excluded. The constraints obtained from 
pair production at LEP and at CDF and from single production at OPAL are also shown. 

Assuming that only one Yukawa coupling h e i is non-vanishingly small, these constraints are 
translated into mass dependent upper limits on the coupling hei, as shown in Fig. 3b-d. 

If the doubly-charged Higgs boson couples only to an electron pair (Fig. 3b) the ee analysis 
rules out H ±± masses below 138 GeV for a coupling h ee of the electromagnetic strength, h ee = 
0.3. The result is compared to the bounds obtained from searches for H ±± pair production at 
LEP [33] and by the CDF experiment [34], and to both the indirect and direct limits obtained by 
the OPAL experiment [9], the latter being the most stringent. The OPAL experiment has also set 
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similar stringent constraints on h ee independently of the Higgs decay mode. These constraints 
also exclude a sizeable H ±± production at HERA via h ee followed by the H ±± decay via /i MM 
or h TT , which is consistent with the non-observation of a resonance signal in the fi/i [31] and 
tt [29] final states in the present HI data. 

Assuming that the doubly-charged Higgs boson couples only to an electron-muon (electron- 
tau) pair, the e/i (er) analysis allows masses below 141 GeV (112 GeV) to be ruled out for 
h e n = 0.3 (h eT = 0.3), as shown in Fig. 3c (Fig. 3d). The HI limits extend the excluded 
region in the electron-muon and electron-tau channels to masses that are beyond those reached 
in previous searches for pair production at LEP [33] and the Tevatron [34]. 



7 Conclusion 

A search for the single production of doubly-charged Higgs bosons coupling to ee, e\x or er 
is presented. In a previous model independent multi-electron analysis, HI observed six events 
with a di-electron mass above 100 GeV, a region where the Standard Model expectation is 
small. Out of the six events, only one is compatible with the signature of a doubly-charged 
Higgs boson. No electron-muon or electron-tau event is found in this mass domain. 

This analysis places new limits on the H ±± mass and its Yukawa couplings h e i to an electron- 
lepton pair. Assuming that the doubly-charged Higgs boson only couples to electron-muon 
(electron-tau) pairs, a limit of 141 GeV (112 GeV) is obtained on the Higgs mass, for a coupling 
h e n = 0.3 (h eT = 0.3) corresponding to an interaction of electromagnetic strength. 
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